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Abstract 

The Bianchi type I cosmological model is brought into a form where the evolution 
of observables is governed by the unconstrained Hamiltonian that coincides with the 
Hamiltonian describing the relative motion of particles in the integrable three-body 
hyperbolic Euler-Calogero-Sutherland system. 
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1 Introduction 

Dealing with cosmological studies, the conceptual problem of the identification of true 
degrees of freedom of gravity, the "golden fleece" of canonical gravity J]J , transforms into 
a real practical problem. In the framework of the conventional formulation of Hamiltonian 
dynamics of the cosmological models (see, e.g., [0, || ||, ||), based on the Dirac || and 
the Arnowitt-Deser-Misner (ADM) treatment of general relativity [0], the solution of the 
problem involves two ingredients: choosing a certain coordinate fixing condition (gauge) 
and solving the constraints. In the application to the minisuperspace models of Bianchi 
type, the standard choice of the intrinsic time gauge leads to a reduced Hamiltonian which 
is in general time- dependent and represents a square- root expression || g, f| ||. In spite 
of the fact that a "relativistic" physicist takes delight in this square-root expression as a 
reminder of the relativistic particle energy, the obtained result is not quite satisfactory. 



1 



Particularly, passing to the quantum theory, it brings the unpleasant difficulties of defining 
the quantized square-root Hamiltonian as an operator in Hilbert space. 

In the present paper, we would like to draw attention to the fact that for the simplest 
Bianchi type I cosmological model, it is possible to use an alternative more attractive 
generator of evolution, free from the above-mentioned difficulties. Below it will be shown 
that one can arrive at a reduced, quadratic in momenta, time-independent Hamiltonian, in 
contrast to the above-mentioned textbook result. Moreover, a transparent correspondence 
between the resulting deparametrized theory and a well-known nonrelativistic many-body 
integrable system naturally arises. Namely, it turns out that the expression for the gen- 
erator of evolution of the unconstrained Bianchi type I cosmology coincides with the 
Hamiltonian describing the relative motion of the particles in the three-body hyperbolic 
Euler-Calogero- Sutherland model. [] 

The plan of the paper is as follows. The first part offers a short review of Bianchi 
cosmological models including as well geometrical aspects as its Dirac generalized dynam- 
ics 



pi] , p4 [23| , [24J]. Then, we restrict our consideration to the simplest Bianchi type I 
cosmology. Analyzing the energy constraint and the extended Poincare-Cartan 1-form, 
we show how to deparametrize the theory, identify the physical degrees of freedom, and 
as a result establish the correspondence with the three-body Euler-Calogero-Sutherland 
model. Finally, using the relation between the dynamics of the Euler-Calogero-Sutherland 
system and the geodesic motion on the space of symmetric matrices, we discuss the general 
solution of the Bianchi type I model. 



2 Bianchi cosmology: geometry and dynamics 

The canonical formulation of the cosmological models intensively uses the geometrical 
symmetry to restrict the gravitational configuration space, the space of all possible pseudo- 
Riemannian metrics defined on a given manifold. In the case of large symmetry of the 
space-time manifold, the gravitational degrees of freedom are effectively reduced to finite 
number and this circumstance essentially relieves the analysis of the theory. Below we 
describe very briefly this effective gravitational configuration space for the homogeneous 
Bianchi type models. For details, we refer to the comprehensive reviews [|. |J . 



2.1 Geometrical settings 

The conventional Hamiltonian analysis requires the existence of a preferred timelike vari- 
able and therefore views the Universe in terms of space plus time. Below in quick review 
of the 3 + 1 geometry of gravitation and its canonical formulation, we follow the definitions 
and notations from the very transparent review by Isenberg and Nester [II]. 



1 This system is a generalization of the many-body Calogero-Sutherland model || [T^, [Tl]; |l2| (see 
also the useful reviews [l4, 15 ), describing particles on a line and interacting via pairwise hyperbolic 
potential l/sinh 2 x, by introduction of spin variables of the particles 16, Il7|. For the other types of 
generalizations by inclusion of spinlike variables in the Calogero-Sutherland model, see, e.g., |l8|, [l9|, p(| . 
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So, the space-time is supposed to be a smooth manifold, {M. , g) endowed with a metric 
g of signature (— , +, +, +), metric-compatible connection, and foliated by a family of one- 
parameter three-dimensional nonintersecting surfaces E t on which the induced metric has 
three positive eigenvalues, i.e., M. = 1 x E f . According to such a space-time foliation, it 
is useful to choose on M. a surface-compatible frame of vector fields (e_i_ , &i, e2, 63) with 
timelike unit-length vector field e_i_ orthogonal to E 4 and three spacelike vector fields f] 
e a — (ei, &2i e3)tangent to it. In the corresponding frame of forms {0 L ,9 a ), the metric g 
reads 

S =-0 ± ®0 ± + 'y ab ff l ®O b , (1) 

with a spatial metric 7 induced on S^. According to this representation, the Lie derivative 
Ce ± , derivative with respect to the proper time along the normal to the hypersurface S t , 
describes the evolution of the dynamical variables. Note that the frame (e_|_,e a ) is not a 
coordinate frame and usually, instead of dealing with Eq. (|I|), the so-called ADM metric 
is exploited. The latter is based on using the coordinate vector fields e = d/dt and 
e a = d/dx a , where the vector field eo is expressed in terms of the normal vector field e_i_ 
and spatial vector field N a e a tangent to the hypersurface £ f as 

e = Ne ± + N a e a . (2) 

However, when investigating the homogeneous Bianchi cosmological models, it is more 
convenient to deal with a special noncoordinate frame instead of the ADM coordinate 
frame for the 3-space. Its construction is dictated by geometrical considerations. By 
definition, in a spatially homogeneous space-time, a three-dimensional Lie group G3 acts 
on the space-time as a group of isometries, such that each orbit on which G 3 acts simply 
transitively is a spacelike hypersurface. The advantage of considering the simply transitive 
action is that we can put the element of G3 into one-to-one correspondence with the 
points of £f and thus the space-time is considered topologically as the product space 
A4 — R x G3. Based on this observation, it is clear that instead of the usual coordinate 
frame of the spatial vector fields we need to choose a new space basis e a , adapted to the 
Lie group structure of the three-dimensional hypersurface S 4 . Namely, the algebra of the 
infinitesimal generators of isometries, i.e., the Killing vector fields £ a , 

C c a Uc (3) 

dictates what kind of basis should be chosen. Indeed, the vector fields £ a provide a basis 
(e , e a ) invariant under the isometries 

£ 5a e = 0, %e a = 0. (4) 

In this case, one can specify the form of the three-dimensional part of the Bianchi metric 

as 

7 = lab u a ®u\ (5) 

2 Following Q, the boldface notation is used to distinguish four-dimensional quantities from three- 
dimensional ones. 

3 Hereafter, the latin indices run over 1,2,3. 
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where uj a are group-invariant 1-forms dual to the vector fields e a with structure coefficients 
C\ c = 2du a (eb, e c ), defined by the structure constants of the homogeneity group C a bc . 

Thus, instead of the ADM coordinate frame, the adapted representation for the 
Bianchi-type metrics in noncoordinate basis is 

g = -(N 2 - N a N a )dt ® dt + 2N a dt ® uj a + ^ ab uj a ® uj b . (6) 

The preferable role of this choice for the frame is transparent: from the Killing equation 

A.S = (7) 

it follows that the functions N, N a , and 7^, entering in the expression for the Bianchi 
metric (^j), depend only on the time parameter t. 



2.2 Generalized dynamics 

As a result of this geometrical analysis, in an adapted basis, all gravitational field config- 
uration variables, the lapse function N, the shift vector N a , and the spatial metric 7^ can 
be identified with a set of ten Lagrangian coordinates of a certain "mechanical system." 
By definition, the dynamics of this system is determined from the Hilbert action 

A= f crR, (8) 

where the ansatz (Q) is plugged into the expression for the space-time scalar curvature R 
and for the four-dimensional volume element cr = yj— g u>° Auo 1 Auo 2 Aw 3 . As a result the 
restricted variational problem (for Bianchi class A models) [] [0, [|] is 



(9) 



*2 ta 

A[N, N a , lab , i ah ] = J dt^N^R - K a a K b b + K ab K ab ) + 2 J ' dt ± (^ a a ) , 

where 3 R is the scalar curvature formed from the spatial metric 7, 

" R =-\ l ah C c da C d ch - \ 7 a y%,C" ac C\ d , (10) 

and 

K ah = ~ [ (7a. C\ c + lbd C d J N C + 7afe ] (11) 

is the extrinsic curvature of the slice E t defined as K ab — — | C e ^ ab . 

The Lagrangian (^) belongs to the class of so-called degenerate ones; its Hessian is 
zero. To deal with its Hamiltonian description, we shall follow the Dirac generalization 
of Hamiltonian dynamics pi , [22] , |23| , p4| . 



4 Bianchi class A models are denned writing the structure constants of the isometry Lie group as 
C% b = ei ab S ld + A[ d 5^ and supposing that C^ d = A a — 0. Here S ab is a three-dimensional second-rank 
symmetric tensor and A a is a three-dimensional vector. 
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Implementing the Legendre transformation on variables N, N a } and ^y ab , we get the 
canonical Hamiltonian 

H C = NH + N a H a (12) 



the primary 

and secondary constraints 

H 



P = 0, P a = (13) 



VT H a b ~ -7T a 7T b 



Vi 3 R = o, (14) 



1 

H a = 2C d ab 7T bc lcd = 0. (15) 



Here the canonical variables (N ,P), (N a ,P a ), and ( , y ab ,n ab ) obey the following nonvan- 
ishing fundamental Poisson brackets relations: 

{N,P} = 1, (16) 
{N a ,P b } = S a b , (17) 

{la b X d } = \{5 c a 5 d b +5 d Jl). (18) 

Due to the reparametrization symmetry of Eq. ([|) inherited from the diffeomorphism 
invariance of the initial Hilbert action, the evolution of the system is unambiguous and it 
is governed by the total Hamiltonian 

H T = NH + N a H a + uP + u a P a , (19) 

with four arbitrary functions u(t) and u a (t). One can verify that the secondary constraints 
are first class and obey the algebra 

{H,H b } = 0, (20) 
{H a , H b } = -C d ab H d ■ (21) 

Thus the dynamics of system @ in the extended phase space, spanned by the canonical 
variables (N , P) , (N a , P a ) , and {'jab ,' n ' ab ), is described by the action rewritten in the 
Hamiltonian form as 

A = f 9 + / d(7c ab lab ) , (22) 



where G is the Poincare-Cartan 1-form 

6 = Ti ab d^/ab + P dN + P a dN a — Ht dt . (23) 
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3 Deparametrized version of the model 



3.1 Misner-Ryan decomposition 

Let us now specialize to the Bianchi type I model. In this case, the three-dimensional 
group of isometries is an Abelian one, i.e., C ab = 0, and thus the momentum constraints 
(115]) are identically satisfied, while the energy constraint reduces to the simplest form 



n = n a \ ab - l -^y h . (24) 

To find the corresponding unconstrained Hamiltonian system with a certain observable 
time parameter, we shall use the variables introduced by Misner || and modified by Ryan 
[J2J. In Misner's representation, the spatial metric 7 is given by 

7« = r 2 e-^% , (25) 

where /3y is a 3 x 3 symmetric traceless matrix and Q is a scalar, both being functions 
of the time parameter only, and r is a constant. f] For our purposes it is convenient, 
following H f|] , to pass to the main-axes decomposition for the nondegenerate symmetric 
matrix /3 and finally write the metric 7 in the Misner-Ryan form as 

e P++V3/3- Q 

- /?'(V) I e^"^- ) /?(\). (20) 
f 

Here R{x) is an orthogonal SO(3, M) matrix parametrized with the three Euler angles Xi- 
From the Jacobian of the transformation (|26|), 

J h*b(tt,J3,x) \ ^ e -m J sinh ( 2 V3/3_) sinh(3/5+ + V3f3.) sinh(3/3+ - y/30_) \ , (27) 




it follows that Eq. (|2q) can be used as a definition of the new configuration variables, 
(Q, 0±), and the three angles (xi, X2, Xs) oia ^Y ^ an eigenvalues of the matrix 7 are different. 
To have a uniqueness of the inverse transformation, we are forced to treat only this type 
of configurations, the so-called principle orbits of the action of the 5*0(3, R) group, while 
the analysis of the orbits with coinciding eigenvalues of the matrix S (singular orbits) 
requires special consideration. 

The point transformation ( p6|) induces the canonical transformation from the coordi- 
nates (j a b , 7i ab ) to the six new canonical pairs (fl ,pn), (fi± ,p±), and (xi ,Pxj)> 

{fi,pn}=l, {/3±,P±} = 1, {Xi,P Xi } = 8ii- (28) 



5 In Eq. (J2^), the constant ro was introduced for convenience in choosing units. For notational 
convenience, in this paper the constant ro is scaled to unity hereafter. 
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Using the requirement of invariance of the symplectic 1-form, 



TT ab d-f ab = Pn dn + P+ d(3+ + p-dfi- + Pxi d Pxi » ( 29 ) 

i=i 

one can find the explicit formulas for the new variables as functions of the old ones and 
finally rewrite the energy constraint as 

2H=^{-pI+pI+pI) 

+1 ( £ + ^ + \ (go) 

4 Vsinh 2 (3/5+ - y/3/3-) sinh 2 (3/5+ + y/3/3-) sinh 2 (2^3/3-)/ 

Note that in this expression all angular variables are gathered into the three right-invariant 
SO(3,R) Killing vector fields r] h 

+ I I Sill Xl /on 

7]i = - sin xi cot X2 P X i + cos Xi P X 2 + ~ Pxs > ( 31 ) 

sm X2 
cos^x 

^2 = cos xi cot x 2 P X i + sm Xi Px2 - Pxs i ( 32 ) 

smx2 

Vs = Pxi , (33) 
satisfying the 50(3, R) Poisson brackets algebra 

{Vh Vj} = UjkVk- (34) 

3.2 Deparametrization and correspondence to the hyperbolic 
three-particle Euler-Calogero- Sutherland model 

In this section, it will be shown that the unconstrained Hamiltonian, describing the evo- 
lution of observables in Bianchi type I cosmology, can be identified with the Hamiltonian 
of the relative motion in the three-particle hyperbolic Euler-Calogero-Sutherland system. 

The direct method to "deparametrize" the initial reparametrization-invariant theory 
with first-class constraints consists of two steps, namely fixation of the gauge and pro- 
jection of the initial degenerate action to the constraint shell. Then the expression for 
the reduced, unconstrained Hamiltonian can be extracted from the Poincare-Cartan 1- 
form. Here it is in order to make one general comment on the gauge-fixing when some 
preferred set of variables is chosen. It turns out that in this case the action projected 
on the constraint shell becomes independent of the part of the variables without impos- 
ing some gauge-fixing condition. These eliminated variables should naturally be treated 
as unphysical ones. Such a phenomenon is well-known from the theory of Hamiltonian 
systems possessing a Lie group symmetry; when the cordinates adapted to the symme- 
try group action are found, then the so-called cyclic (or ignorable) coordinates disappear 
from the effective Hamiltonian. In the reparametrization-invariant theories, there is a 
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certain peculiarity of this effect. Namely, when such a set of well-defined variables exists, 
then after reduction, one of these cyclic variables takes the role of time in the corre- 
sponding deparametrized system. The existence of such an exceptional variable in the 
reparametrization-invariant theories means that the equivalent unconstrained system can 
be written in an autonomous form. 

Below we shall see that the model under consideration can be brought into a class of 
such Hamiltonian systems. [] To see how the proper variable which takes the role of time 
in the unconstrained system appears, let us consider the energy constraint (^) in terms 
of the Misner-Ryan coordinates. One can easily verify that performing the canonical 
transformation 

T = 6 p\> Ut = Y2 P ^ (35) 
the energy constraint (|30"D linearizes with respect to the momentum LPr, [] 

2n = -U T + h((3 ± ,p ± , Vi ), (36) 

where h(/3± ,p± , rji) is given by 

i 1 ( | ^ | A (37) 

4 Vsinh 2 (3/3+ - y/3/3-) sinh 2 (3/3+ + v 7 ^-) sinh 2 (2 v / 3/3-) / 

Note that after all of these transformations, the Poincare-Cartan 1-form (^) changes by 
a total differential, 

3 1 
6 = U T dT + P+ d(3 + + p_ dp. + P^ d Xi + PdN + P a dN a - H T dt - - d(P n Q) . (38) 

i=i 

Now to find the evolution parameter and the Hamiltonian of the deparametrized the- 
ory, we project Eq. flSlf ) to the constraint shell P = P a = H = 7i a = 0, Q 

3 

\constraints = P+ d P+ + V-^- + ^Pxi^Xi + K(3 ± , P± , Vi )dT . (39) 

i=l 

The expression for the projected 1-form (|39|) prompts the following interpretation. Sup- 
posing that the deparametrized version of the Bianchi type I model has a conventional 
Poincare-Cartan 1-form, 

3 

Q* BI = p+d(3+ + p-d/3- + Y Px* d Xi ~ H B i(J3± , p± , rj^dr , (40) 

i=l 

6 The corresponding representation for the Friedmann cosmological model with scalar and spinor fields 
has been discussed recently in p6[ . 

7 Such a linearity in one of the momenta forms of the constraints in accordance with Kuchaf 's obser- 
vation |^5| shows the possibility of a global deparametrization of the degenerate theory. 

8 We omit all differentials assuming that proper boundary conditions are taken to derive the classical 
equation of motion. 
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we make the following identifications: the initial dynamical variable T(t) "metamor- 
phosis" to the evolution parameter r = — T of the corresponding deparametrized sys- 
tem and the Hamiltonian, as a time-independent generator of evolution, is given by 
Hbi = h((3± , p± , T]i) . Thus to summarize, we arrive at the unconstrained Hamiltonian 
for the Bianchi type I model in the form 

H BI = ^(pl+p 2 _) 

1 / 7/2 1)1 Vl 



+ 1 Jl + ^ + (41) 

4 Vsinh 2 (3/3 + - V3[3„) sinh 2 (3/3+ + V3/3-) sinh 2 (2 v / 3/?-)/ 

One final comment is in order concerning the form of the derived Hamiltonian. As 
was mentioned in the Introduction, the appearance of the square-root expression for the 
unconstraned Hamiltonian of the Bianchi cosmologies is usually interpreted in a certain 
analogy to relativistic particles 0, |3J]. In contrast to this picture, our analysis points to 
the existence of a certain correspondence with a system of nonrelativistic particles. 

To describe this correspondence, let us recall the formulation of the integrable three- 
particle hyperbolic Euler-Calogero-Sutherland (ECS) model (16|, |I7|| . The particles are 
characterized by two types of variables. The first ones are the canonical coordinates 
(Pi ,pi) describing the position and momenta, 

{Pi,p j } = 6 ij , (42) 

and the other are the so-called "internal" coordinates l a b = —ha, obeying the SO(3, R) 
Poisson brackets algebra, 

{lab, ha 1 } = SaJbd ~ ^adhc + <Wac ~ ^bJad ■ (43) 

The dynamics of both variables is determined by the Hamiltonian 

flra = £rf + 5E_4_W- (44) 



(»#) 



Now it is easy to find the relation between the derived unconstrained Hamiltonian 
y]) and the Euler-Calogero-Sutherland Hamiltonian (f44|). To achieve this, we identify 
the internal variables with the Killing vectors rji of the Bianchi type I model, 

kj = tijkVk , (45) 

while the variables (j3±, p±) are identified with the relative variables in the Jacobi system 
for (fa, pi), 

/3 1 = P + + V3P..-X, Pl = lp + + -^-=p^-\p, (46) 

6 2V3 3 

P 2 = P + -V3P--X, p2 = ±p + --L p _-}-P, (47) 

P 3 = -2P + -X, p 3 = _I p+ _Ip. (48) 
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Here (X, P) are the corresponding center-of-mass coordinates. Now one can verify that 
the expression (f44|) takes the following separable form: 

Hecs = \p 2 + Hbi (49) 
o 

and thus one arrives at the above-stated correspondence between the Hamiltonian, de- 
scribing the relative motion in a three-particle hyperbolic Euler-Calogero-Sutherland sys- 
tem, and the deparametrized version of the Bianchi type I model Hamiltonian Hbi- 



3.3 Solution of the unconstrained equations of motion 

The established relation to the many-body integrable model allows us to write down the 
general solution for the Bianchi type I metric. Here it is in order to note that the solution 
of the initial Einstein equations for the Bianchi-type metrics contain arbitrary functions, 
reflecting the reparametrization invariance, of the theory and the solution should be clas- 
sified according to their equivalence. Passing to the reduced system, when the constraints 
generating this reparametrization are eliminated, this freedom is fixed. However, if the 
system has the additional rigid symmetry, not all of solutions describe the nonequivalent 
metrics and now the problem of classification reduces to the identification of the metrics, 
with a fixed number of parameters. Below we shall present the general solution for the 
reduced Bianchi type I model, while the detailed analysis of the essential parameters and 
their origin will be done elsewhere. 

To write down an explicit form of the solution for the general metric of Bianchi type 
I space-time, we recall the roots of the integrability of the Euler-Calogero-Sutherland 
model. It is known (see, e.g., [19|, [HJ and references therein) that the integral curves of 



the Hamiltonian ([44]) are in one-to-one correspondence with the geodesies of the space of 
symmetric matrices S endowed with the bi-invariant metric, 



J ECS 



~ J TriS-US) 2 . (50) 



To prove this dynamical equivalence, it is enough to use again the main-axes decomposi- 
tion for the symmetric matrix S of the form 

S = R T e 2X R (51) 

with the diagonal matrix X = diag||/5i, /? 2 , /5 3 ||. Assuming that SV,- are six independent 
Lagrange coordinates, from the action ( |50| ) it follows that the equation of motion can be 
written in the form of conservation law, 

" ? =0, (52) 



dt \ dt 

and therefore the general solution for the matrix S can be represented in a simple form, 

S(t) = S e Jt . (53) 
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Here Sq and J are arbitrary constant matrices, Sq is a symmetric matrix, the value of the 
matrix S(t) at an initial moment t = 0, while J is an integral of motion, J = Sq Sq. [] 
Comparison of the two expressions (|51| ) and ([53]) gives the solution of the Euler-Calogero- 
Sutherland model and therefore of the Bianchi type I metric due to the correspondence 
derived in the previous section. 



4 Concluding remarks 

Following the Dirac generalized Hamiltonian approach we derive the deparametrized ver- 
sion of the Bianchi type I cosmology in the form of the integrable three-particle Euler- 
Calogero-Sutherland model. It is worth posing the question of the possibility to extend the 
result presented here and to find the quadratic time-independent unconstrained Hamil- 
tonians for the other more complicated Bianchi systems. In seeking the answer, first of 
all it is necessary to investigate the rigid symmetries that Bianchi models possess, be- 
cause the additional rigid symmetry as a rule leads to the conception of preferred time 
choice |[25|| . More precisely, the existence of a timelike Killing vector field allows us to 
define an intrinsic time and construct the observables, whose dynamics is governed by a 
time-independent Hamiltonian. 

Another application of the discussed correspondence is to analyze the metric of Bianchi 
type I cosmology with respect to its dependence on the essential parameters. In a forth- 
coming publication, we plan to discuss in detail a certain parametrization of the metric 
using the solution flo"3"| ) of the three-body Euler-Calogero-Sutherland problem and we shall 
attempt to clarify how the topology of the flat Bianchi type I cosmology emerges. One 
can find the recent discussions of the topological characteristics of the three-dimensional 
space in the cosmology in review pTI . 
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